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Abstract 



Branching fractions of the decays Hi, — )• HcTt 7r"'"7r relative to Hi, — )• Heir are 
presented, where Hb {H^) represents ^° (£»+), B' (L>°), B^ (L>+) and A° (A+). 
The measurements are performed with the LHCb detector using 35 pb~^ of data 
collected at ^/s = 7 TeV. The ratios of branching fractions are measured to be 



= 2.38 ±0.11 ±0.21 
= 1.27 ±0.06 ±0.11 
-- 2.01 ±0.37 ±0.20 
1.43 ±0.16 ±0.13. 
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We also report measurements of partial decay rates of these decays to excited charm 
hadrons. These results are of comparable or higher precision than existing measure- 
ments. 
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1 Introduction 



Over the last two decades, a wealth of information has been accumulated on the decays 
of 6-hadrons. Measurements of their decays have been used to test the CKM mecha- 
nism [T] for describing weak decay phenomena in the Standard Model, as well as provide 
measurements against which various theoretical approaches, such as HQET [2] and the 
factorization hypothesis, can be compared. While many decays have been measured, a 
large number remain either unobserved or poorly measured, most notably in the decays 
of Bg mesons and A° baryons. Among the largest hadronic branching fractions are the 
decays Hb H^tt-tt+tt-, where Ht (H^) represents 5° B~ 5° {D+) and 

(A+). The first three branching fractions were determined with only 30-40% accuracy, 
and the A° — A^tt^tt+tt^ branching fraction was unmeasured. 

Beyond improving our overall understanding of hadronic b decays, these decays are 
of interest because of their potential use in CP violation studies. It is well known that 
the Cabibbo-suppressed decays B' DK' [31 H E] and B° -)> DfK^ [7] provide 
clean measurements of the weak phase 7 through time-independent and time-dependent 
rate measurements, respectively. Additional sensitivity can be obtained by using B^ — )■ 
D'^Tx~ [8] decays. As well as these modes, one can exploit higher multiplicity decays, 
such a.s_W DK*°, B~ DK~tx+ti- [9] and B° ^ DfK^n^n^. Moreover, the 
decay 5° — ?■ D^tc^ti^ti" has been used to measure Am^ [ID], and with a sufficiently large 
sample, provides a calibration for the flavor-mistag rate for the time-dependent analysis 
of Bl DfK^ir^TT^. 

The first step towards exploiting these multi-body decays is to observe them and 
quantify their branching fractions. The more interesting Cabibbo-suppressed decays are 
O(A^) in the Wolfenstein parameterization [TT], and therefore require larger data samples. 
Here, we present measurements of the Cabibbo-favored Hfy — )■ Hcn~7T~^TT~ decays. The 
leading amplitudes contributing to these final states are shown in Fig. [H Additional 
contributions from annihilation and VT-exchange diagrams are suppressed and are not 
shown here. Note that for the B~ and A^ decays, unlike the 5° and 5^, there is potential 
for interference between diagrams with similar magnitudes. In Ref . [12] , it is argued that 
this interference can explain the larger rate for B~ — j- D^tt~ compared to B^ — )■ D~^7r~ . 
Thus, it is interesting to see whether this is also true when the final state contains three 
pions. 

In this paper, we report measurements of the Hb — )■ HcTT~tt^it~ branching frac- 
tions, relative to — )■ Hcn~. We also report on the partial branching fractions, 
Hb H*TT~, H* HcTT^n", where Hb is either 5°, B~, or A^, and H* refers to 
Di(2420)+'°, £'^(2460)°, Ac(2595)+, or A,(2625)+. We also present results on the partial 
rates for A^ — Sc(2544)°'+"'~7r^7r~. Charge conjugate final states are implied throughout. 

2 Detector and Trigger 

The data used for this analysis were collected by the LHCb experiment during the 2010 
data taking period and comprise about 35 pb~^ of integrated luminosity. LHCb has 
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Figure 1: Feynman diagrams for iJ;, — Hcji" and if^ — )■ Hc7i~n~^7i~ decays. Figs, (a) and 
(b) show external tree diagrams, (c) and (d) show color-suppressed tree diagrams {B~ and 
A° only), and (e) shows the Cabibbo-suppressed external tree diagram, only accessible to 
meson. 

excellent capabilities to trigger on and reconstruct bottom and charm hadrons. The most 
important element of the detector for this analysis is a charged particle tracking system 
that covers the forward angular region from about 15 — 350 mrad and 15 — 250 mrad in 
the horizontal and vertical directions, respectively. It includes a 21 station, one- meter 
long array of silicon strip detectors (VELO) that come within 8 mm of the LHC beams, 
a 4 Tm dipole magnetic field, followed by three multi-layer tracking stations (T-stations) 
downstream of the dipole magnet. Each T-station is composed of a four layer silicon 
strip detector (IT) in the high occupancy region near the beam pipe, an eight layer straw 
tube drift chamber (OT) composed of 5 mm diameter straws outside this high occupancy 
region. Just upstream of the dipole magnet is a four-layer silicon strip detector (TT). 
Overall, the tracking system provides an impact parameter (IP) resolution of ~ 16/im + 
SOfim/pT (transverse momentum, pt in GeV/c), and a momentum resolution that ranges 
from ap/p ~ 0.4% at 3 GeV/c to ~ 0.6% at 100 GeV/c. Two Ring Imaging Cherenkov 
Counters (RICH) provide a kaon identification efficiency of ~95% for a pion fake rate of a 
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few percent, integrated over the momentum range from 3—100 GeV/c. Downstream of the 
second RICH is a Preshower/Scintillating Pad Detector (PS/SPD), and electromagnetic 
(ECAL) and hadronic (HCAL) calorimeters. Information from the ECAL/HCAL is used 
to form the hadronic triggers. Finally, a muon system consisting of five stations is used 
for triggering on and identifying muons. 

To reduce the 40 MHz crossing rate to about 2 kHz for permanent storage, LHCb 
uses a two-level trigger system. The first level of the trigger, LO, is hardware based and 
searches for either a large transverse energy cluster [Et > 3.6 GeV) in the calorimeters, or 
a single high pt or di-muon pair in the muon stations. Events passing LO are read out and 
sent to a large computing farm, where they are analyzed using a software-based trigger. 
The first level of the software trigger, called HLTl, uses a simplified version of the offline 
software to apply tighter selections on charged particles based on their px and minimal IP 
to any primary vertex (PV), defined as the location of the reconstructed pp collision(s). 
The HLTl trigger relevant for this analysis searches for a single track with IP larger than 
125 /im, Pt > 1.8 GeV/c, p > 12.5 GeV/c, along with other track quality requirements. 
Events that pass HLTl are analyzed by a second software level, HLT2, where the event is 
searched for 2, 3, or 4-particle vertices that are consistent with fe-hadron decays. Tracks 
are required to have p > 5 GeV/c, pt > 0.5 GeV/c and IP larger than 16 to any PV, 
where the value is obtained assuming the IP is equal to zero. We also demand that 
at least one track has px > 1.5 GeV/c, a scalar p-r sum of the track in the vertex exceed 
4 GeV/c, and that the corrected masqj is between 4 and 7 GeV/c^. These HLT trigger 
selections each have an efficiency in the range of 80—90% for events that pass typical 
offline selections for a large range of B decays. A more detailed description of the LHCb 
detector can be found in Ref. [T3] . 

Events with large occupancy are known to have intrinsically high backgrounds and to 
be slow to reconstruct. Therefore such events were suppressed by applying global event 
cuts (GECs) to hadronically triggered decays. These GECs included a maximum of 3000 
VELO clusters, 3000 IT hits, and 10,000 OT hits. In addition, hadron triggers were 
required to have less than 900 or 450 hits in the SPD, depending on the specific trigger 
setting. 

3 Candidate Reconstruction and Selection 

Charged particles likely to come from a 6-hadron decay are first identified by requiring 
that they have a minimum IP with respect to any PV of more than 9. We also 
require a minimum transverse momentum, px > 300 MeV/c, except for Hf, — )■ iJc7r^7r+7r~ 
decays, where we allow (at most) one track to have 200 < px < 300 MeV/c. Hadrons are 
identified using RICH information by requiring the difference in log-likelihoods (ALL) 
of the different mass hypotheses to satisfy ALL{K — vr) > —5, ALL(p — vr) > —5 and 

^The corrected mass is defined as Mcor = \/ + Ptrans > where M is the invariant mass of the 2, 3 
or 4-track candidate (assuming the kaon mass for each particle), and ptrans is the momentum imbalance 
transverse to the direction of flight, deflned by the vector that joins the primary and secondary vertices. 
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ALL{K — tt) < 12, for kaons, protons and pions, respectively. These particle hypotheses 
are not mutually exclusive, however the same track cannot enter more than once in the 
same decay chain. 

Charm particle candidates are reconstructed in the decay modes D'^ — )■ K~7r~^, — t- 
K~7r~^n~^, — )■ K~^K~7r~^ and — >■ pK~7r~^. The candidate is associated to one of the 
PVs in the event based on the smallest IP between the charm particle's reconstructed 
trajectory and all PVs in the event. A number of selection criteria are imposed to reduce 
backgrounds from both prompt charm with random tracks as well as purely combinatorial 
background. To reduce the latter, we demand that each candidate is well separated from 
the associated PV by requiring that its flight distance (FD) projected onto the z-axis is 
larger than 2 mm, the FD > 4dfl, and that the distance in the transverse direction (Ai?) 
is larger than 100 /im. Background from random track combinations is also suppressed 
by requiring the vertex fit x^/ndf< 8, and pr > 1-25 GeV/c (1.5 GeV/c for in 

— Dfn^.) To reduce the contribution from prompt charm, we require that the 
charm particle has a minimal IP larger than 80 fim and IP > 12.25 with respect 
to its associated PV. For — )■ K~^K~7r~^, we employ tighter particle identification 
requirements on the kaons, namely ALL{K — vr) > 0, if the K^K~ invariant mass is 
outside a window of ±20 MeV/c^ of the (f) mass [15]. Lastly, we require the reconstructed 
charm particles masses to be within 25 MeV/c^ of their known values. 

The bachelor pion for Hi, — ?■ Hc7r~ is required to have pt > 0.5 GeV/c, p > 5.0 GeV/c 
and IP > 16. For the Svr vertex associated with the Hf, — > HcH~tt~^7i~ decays, we apply 
a selection identical to that for the charm particle candidates, except we only require the 
Pt of the Sir system to be larger than 1 GeV/c and that the invariant mass is in the range 
from 0.8 GeV/c^ < M{m) < 3.0 GeV/c^. 

Beauty hadrons are formed by combining a charm particle with either a single pion 
candidate (for Hh — )■ Hc7r~) or a Svr candidate (for Hj, — >■ Hcir^TT^Tr' .) The 6-hadron is 
required to have a transverse momentum of at least 1 GeV/c. As with the charm hadron, 
we require it is well-separated from its associated PV, with FD larger than 2 mm, FD 
x"^ > 49 and Ai? > 100 fim. We also make a series of requirements that ensure that the 
6-hadron candidate is consistent with a particle produced in a proton-proton interaction. 
We require the candidate to have IP<90 fim, IP < 16, and that the angle 6 between 
the 6-hadron momentum and the vector formed by joining the associated PV and the 
decay vertex satisfies cos^ > 0.99996. To ensure a good quality vertex fit, we require a 
vertex fit x^/ndf< 6 (8 for Hb HcTx^ .) 

To limit the timing to process high occupancy events, we place requirements on the 
number of tracks^ in an event. For S° D+vr" and S° I^+tt", the maximum number 
of tracks is 180, and for A° — A^vr" and B~ — >■ D^ti~ it is 120. These selections are 99% 
and 95% efficient, respectively, after the GECs. The — >■ Hf{K~Ti^TT~ selection requires 
fewer than 300 tracks, and thus is essentially 100% efficient after the GECs. 

Events are required to pass the triggers described above. This alone does not imply 

^This is the with respect to the FD=0 hypothesis. 

''Here, tracks refer to charged particles that have segments in both the VELO and the T-stations. 
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that the signal fe-hadron decay was directly responsible for the trigger. We therefore also 
require that one or more of the signal 6-hadron daughters is responsible for triggering the 
event. We thus explicitly select events that Triggered On the Signal decay (TOS) at LO, 
HLTl and HLT2. For the measurements of excited charm states, where our yields are 
statistically limited, we also make use of LO-triggers that Triggered Independently of the 
Signal decay (TIS). In this case, the LO trigger is traced to one or more particles other 
than those in the signal decay. 

Lastly, we note that in Hb — )■ iJc7r^7r+7r~ candidate events, between 4% and 10% have 
multiple candidates (mostly two) in the same event. In such cases we choose the candidate 
with the largest transverse momentum. This criterion is estimated to be (75 ± 20)% 
efficient for choosing the correct candidate. For Hf, — )■ Hc7t~ multiple candidates occur in 
less than 1% of events, from which we again choose the one with the largest pt- 

3.1 Selection Efficiencies 

Selection and trigger efficiencies are estimated using Monte Carlo (MC) simulations. The 
MC samples are generated with an average number of interactions per crossing equal 
to 2.5, which is similar to the running conditions for the majority of the 2010 data. 
The fe-hadrons are produced using pythia [16] and decayed using evtgen [17]. The 
Hf, — )■ Hcir'n^TT' decays are produced using a cocktail for the tttttt system that is 
~2/3 ai(1260)" — )■ p°7r~ and about 1/3 non-resonant p^n~. Smaller contributions from 
D?(2420)7r and D*°(2460)7r are each included at the 5% level to B~ D^n'TV+n^ and 
2% each for D^tt^tt^tt^ . For A° — A^tt^tt^tt", we include contributions from 

Ac(2595)+ and Ac(2625)+, which contribute 9% and 7% to the MC sample. The detector 
is simulated with geant4 [18], and the event samples are subsequently analyzed in the 
same way as data. 

We compute the total kinematic efficiency, e^in from the MC simulation as the fraction 
of all events that pass all reconstruction and selection requirements. These selected events 
are then passed through a software emulation of the LO trigger, and the HLT software used 
to select the data, from which we compute the trigger efficiency (ctrig). The efficiencies 
for the decay modes under study are shown in Table [T] Only the relative efficiencies are 
used to obtain the results in this paper. 

4 Reconstructed Signals in Data 

The reconstructed invariant mass distributions are shown in Figs. |2] and [3] for the signal 
and normalization modes, respectively. Unbinned likelihood fits are performed to extract 
the signal yields, where the likelihood functions are given by the sums of signal and several 
background components. The signal and background components are shown in the figures. 
The signal contributions are each described by the sum of two Gaussian shapes with equal 
means. The relative width and fraction of the wider Gaussian shape with respect to the 
narrower one are constrained to the values found from MC simulation based on agreement 
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Table 1: Summary of efficiencies for decay channels under study. Here, ekin is the total 
kinematic selection efficiency, etrig is the trigger efficiency, and etot is their product. The 
uncertainties shown are statistical only. 

Decay ekin etrig etot 
(%) (%) (%) 

E° ^ D+7r"7r+7r~ 0.153 ±0.003 22.6 ± 0.5 0.0347 ± 0.0011 

5" ^ D°7r-7r+7r- 0.275 ± 0.007 27.4 ± 0.6 0.0753 ± 0.0019 

5° ^ D+TT-TT+TT- 0.137 ±0.003 24.9 ±0.7 0.0342 ± 0.0012 

^ A+TT-TT+TT- 0.110 ±0.005 24.0 ±0.7 0.0264 ± 0.0008 

W ^ D+Ti- 0.882 ±0.014 20.8 ±0.3 0.184 ±0.004 

B- ^ DH- 1.54 ±0.02 27.4 ±0.3 0.421 ± 0.007 

B^^^D+TT- 0.868 ±0.010 23.1 ±0.2 0.201 ± 0.003 

A°^A+7r- 0.732 ±0.015 24.7 ±0.4 0.181 ± 0.004 



with data in the large yield signal modes. This constraint is included with a 10-12% 
uncertainty (mode-dependent), which is the level of agreement found between data and 
MC simulation. The absolute width of the narrower Gaussian is a free parameter in the 
ffi, since the data shows a slightly worse (~10%) resolution that MC simulation. 

For B^ — 7- Dfn^ and B^ — )■ D^tt^tt+tt^ decays, there are peaking backgrounds from 
B^ — > D~^7r~ and B^ — )■ D^tt^ti^tt^ just below the B^ mass. We therefore fix their core 
Gaussian widths as well, based on the resolutions found in data for the kinematically 
similar B^ — >■ D^tt" and B^ — )■ D^h'tt^h' decays, scaled by 0.93, which is the ratio of 
expected widths obtained from MC simulation. 

A number of backgrounds contribute to these decays. Below the 6-hadron masses there 
are generally peaking background structures due to partially reconstructed B decays. 
These decays include -B(s) — )■ D*^-|7r(7r7r), with a missed photon, 7r° or tt"*", as well as 

— D(^s^p~ , where the vr" is not included in the decay hypothesis. For the B^ — )■ D'^7t~ 
and B~ — D^vr" decays, the shapes of these backgrounds are taken from dedicated signal 
MC samples. The double-peaked background shape from partially-reconstructed D*7r 
decays is obtained by fitting the background MC sample to the sum of two Gaussian 
shapes with different means. The difference in their means is then fixed, while their 
average is a free parameter in subsequent fits to the data. For B^ — )■ D~^'ir~7r~^7r~ and 
B~ — )■ D°7r~7r"'"7r~, the shape of the partially-reconstructed D*7r7r7r background is not as 
easily derived since the helicity amplitudes are not known. This low mass background is 
also parametrized using a two-Gaussian model, but we let the parameters float in the fit to 
the data. For 5° — )■ -D+vr" and 5° — )■ Dfir^-K'^'iT^, we obtain the background shape from a 
large B^ — D^X inclusive MC sample. Less is known about the A° hadronic decays that 
would contribute background to the A° A^vr" and A° — A+7r~7r+7r~ invariant mass 
spectra. For A° — )■ A+7r~7r"'"7r~, we see no clear structure due to partially- reconstructed 
backgrounds. For A° — > A+tt", there does appear to be structure at about 5430 MeV/c^, 



6 





Mass (MeV/c^) Mass (MeV/c^) 



Figure 2: Invariant mass distributions for — )■ D^tt^tt+tt^ (top left), -B" — ?■ D'^tt^tt^tt" 
(top right), 5° D+TT^TT+TT" (bottom left), and A° A+tt^tt+tt^ (bottom right). Fits 
showing the signal and background components are indicated, and are described in the 
text. 

which may be due to A+p^. The enhancement is described by a single Gaussian above 
the combinatoric background, which, given the limited number of events, provides a good 
description of this background. 

There are also so-called reflection backgrounds, where fully reconstructed signal decays 
from one 6-hadron decay mode produce peaking structures in the invariant mass spectra 
of other decay modes when one of the daughter particles is misidentified. For B — t- 
Dtt^ {tt^tt^), there are reflections from B — t- DK^ (tt^tt^) Cabibbo-suppressed decays, 
where the kaon is misidentified as a pion. Due to the Cabibbo suppression and the 
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Figure 3: Invariant mass distributions for 5° — > D^n^ (top left), — t- D^ti^ (top right), 
5^ — 7- D^-K' (bottom left), and A° — t- A+tt" (bottom right). Fits showing the signal and 
background components are indicated, and are described in the text. 

excellent RICH performance, their contributions are limited to the 1% level. The shape 
of this misidentification background is taken from MC simulation and is constrained to 
be (1 ± 1)% of the signal yield. 

For the 5° — t- Dfn^ and — )■ D^tt^tt+tt" decays, there are reflection backgrounds 
from B^ — )■ D^7T~ and B^ — D+tt^tt+tt^ modes, when either of the tt"*" from the 
decay is misidentified as a K^. This cross- feed background is evaluated in two ways. 
First, we take our B^ — )■ D^n^ {B^ — )■ D^7T~7r~^7r~) data, which have very loose particle 
identification (PID) requirements on the pions, and apply the kaon PID selection to them. 
If either of the two pions pass, and the recomputed {KKir) mass is within the mass 
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window, the candidate is counted as a reflection background. Using this technique, we 
find (5.3 ±0.4)% [(6.3 ±0.6)%] of D+tt^ [5° ^ D+TT~7r+7r~] signal decays refiect 

into the — D^it~ [5^ — Dfn^n^n^] signal region. In the second method, we 
apply a yr-faking-ii' misidentification matrix (in bins of p and px), obtained from a D*~^ 
data calibration sample to the B^ — >■ D^7i~ (or 5° — )■ D^ti~tt^ti~) signal MC sample, 
followed by the mass window requirement (after replacing the pion mass with the 
kaon mass.) The results of this second procedure are (4.4 ± 0.3)% for B^ — D^7r~ and 
(5.2±0.4)% for B^ — t- D+tt^tt+tt", both of which are consistent with the first method. We 
therefore constrain the peaking background from B^ — )■ D~^7t~ [B^ — t- D^ti^tt^ti^) into 
5° ^ D+TT- (5° ^ D+TT-TT+TT-) to be (4.0 ± 1.5)% ((5.0 ± 2.0)%), where the Gaussian 
constraint is conservatively assigned a 40% relative uncertainty. The shape of this peaking 
background is obtained from MC simulation and is well-described by a single Gaussian of 
mean 5350 MeV/c^ and width 30 MeV/c^. This shape is in good agreement with what is 
observed in data. 

The second refiection background to S° — )■ -D^vr" (_B° — )■ D+vr'Tr+vr") is A° — )■ A+vr" 
(A^ — A+7r~7r+7r~), where the proton from the A^ decay is misidentified as a kaon. This 
is similar to the B^ refiection, except here the A° yield is significantly smaller, obviating 
the need for making an explicit ALL{K — p) requirement to reject protons. The A° 
refiection background is evaluated using the first technique as described above leading 
to refiection rates of (15 ± 3)% for A+tt" into Dfn- and (20 ± 4)% for 

A° A+7r~7r"'"7r~ into B^ — )■ D+7r~7r~^7r~. We conservatively assign a 20% uncertainty 
on this rate based on the agreement between data and MC simulation. The asymmetric 
shape of this background is described by the simulation, which is consistent with the 
shape observed in data. The combinatorial background is modeled with an exponential 
distribution. The fits are superimposed on the data in Figs. |2] and [3l and the fitted yields 
are summarized in Table [2J 

The ratios of branching ratios are given by: 



Slg 

tot 



where the Y factors are the observed yields in the signal and normalization modes, and 
etot are the total selection efficiencies. 



5 Systematic Uncertainties 

Several sources contribute uncertainty to the measured ratios of branching fractions. Be- 
cause we are measuring ratios of branching fractions, most, but not all of the potential 
systematics cancel. Here, we discuss only the non- cancelling uncertainties. With regard 
to the reconstruction of the — t- HcT^^-k^'k^ and Hi^ — )■ HcTi^ decays, the former has 
two additional pions which need to pass our selections, and the Sir system needs to pass 
the various vertex-related selection criteria. The track reconstruction efficiency and un- 
certainty are evaluated by measuring the ratio of fully reconstructed J/ip^s to all J/ip^s 
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Table 2: Summary of yields for the branching fraction computation. Uncertainties are 
statistical only. 



Decay 




Yield 


Decay 


Yield 


5° ^ D + TT 


TT + TT 


1150 ±43 




2745 ± 66 


B- D% 




950 ± 41 


B- D^n- 


4244 ± 90 






138 ± 23 




434 ± 32 


AO A+TT- 




174 ± 18 


A" ^ A+TT- 


853 ± 36 



obtained from an inclusive single muon trigger, where only one of the muons is required to 
be reconstructed. After reweighting the efficiencies to match the kinematics of the signal 
tracks, the uncertainty is found to be 2% per track, which leads to a 4% uncertainty in 
the branching fraction ratios. The IP resolution in data is about 20% worse than in the 
simulation, leading to (i) a larger efficiency for tracks to pass the IP-related cuts (as well 
as larger background), and (ii) a lower efficiency to pass the vertex selections, for data 
relative to the value predicted by simulation. The first of these is studied by reducing 
the IP requirement in simulation by 20%, and the second by smearing the vertex 
distribution in simulation until it agrees with data. The combined correction is found to 
be 1.02 ±0.03. 

Another potential source of systematic uncertainty is related to the production and de- 
cay model for producing the Hcimn final state. We have considered that the px spectrum 
of the pions in the Sir system may be different between simulation and data. To estimate 
the uncertainty, we reweight the MC simulation to replicate the momentum spectrum 
of the lowest momentum pion (amongst the pions in the Sn vertex.) We find that the 
total efficiency using the reweighted spectra agrees with the unweighted spectra to within 
3%. We have also investigated the effect of differences in the px spectra of the charm 
particle, and find at most a 1% difference. Our candidate selection is limited to the mass 
region M(7r7r7r) < 3 GeV/c^. Given that the phase space population approaches zero as 
M(7r7r7r) — t- 3.5 GeV/c^ {i.e., Mb — Mo) and that the simulation reasonably reproduces 
the 7r~7r+7r~ mass spectrum, we use the simulation to assess the fraction of the vrvrvr mass 
spectrum beyond 3 GeV/c^. We find the fraction of events above 3 GeV/c^ is (3.5 — 4.5)% 
for the decay modes under study. We apply a correction of 1.04 ± 0.02, where we have 
assigned half the correction as an estimate of the uncertainty. In total, the correction for 
production and decay models is 1.04 ± 0.04. 

As discussed in Sec. [3l we choose only one candidate per event. The efficiency of this 
selection is estimated by comparing the signal yield in multiple-candidate events before 
and after applying the best candidate selection. The selection is estimated to be (75±20)% 
efficient. In the Hh — t- HcTt^ti^h^ the multiple candidate rate varies from 4% to 10%, so 
we have corrections that vary from 1.01 to 1.03. For Hb — > Hcir^ , this effect is negligible. 
The corrections for each mode are given in Table [31 

For the trigger efficiency, we rely on signal MC simulations to emulate the online 
trigger. The stability of the relative trigger efficiency was checked by reweighting the 
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fe-hadron pj- spectra for both the signal and normaUzation modes, and re-evaluating the 
trigger efficiency ratios. We find maximum differences of 2% for LO, 1% for HLTl and 1% 
for HLT2, (2.4% total) which we assign as a systematic uncertainty. 

Fitting systematics are evaluated by varying the background shapes and assumptions 
about the signal parameterization for both the Hj, — )■ ifc7r~vr+7r~ and Hh — )■ Hc7r~ modes 
and re-measuring the yield ratios. For the combinatorial background, using first and 
second order polynomials leads to a 3% uncertainty on the relative yield. Refiection 
background uncertainties are negligible, except for — )■ D+tt^tt+tt" and 5° — )■ Df7i~, 
where we find deviations as large as 5% when varying the central value of the constraints 
on the — >■ D+7r~7r+7r~ and — D+vr" refiections by ±1 standard deviation. We 
have checked our sensitivity to the signal model by varying the constraints on the width 
ratio and core Gaussian area fraction by one standard deviation (2%). We also include 
a systematic uncertainty of 1% for neglecting the small radiative tail in the fit, which is 
estimated by comparing the yields between our double Gaussian signal model and the 
sum of a Gaussian and Crystal Ball [20] line shape. Taken together, we assign a 4% 
uncertainty to the relative yields. For the B^ branching fraction ratio, the total fitting 
uncertainty is 6.4%. 

Another difference between the Hi, — )■ HcH~ and Hb — HcH^-k^-k^ selection is the 
upper limit on the number of tracks. The efficiencies of the lower track multiplicity 
requirements can be evaluated using the samples with higher track multiplicity require- 
ments. Using this technique, we find corrections of 0.95±0.01 for the B~ and A|) branching 
fraction ratios, and 0.99 ± 0.01 for the 5° and 5° branching fraction ratios. 

We have also studied the PID efficiency uncertainty using a D*^ calibration sample 
in data. Since the PID requirements are either common to the signal and normalization 
modes, or in the case of the bachelor pion(s), the selection is very loose, the uncertainty 
is small and we estimate a correction of 1.01 ± 0.01. We have also considered possible 
background from Hi, — )■ H(.D~ which results in a correction of 0.99 ± 0.01. 

All of our MC samples have a comparable number of events, from which we incur 3-4% 
uncertainty in the efficiency ratio determinations. The full set of systematic uncertainties 
and corrections are shown in Table |3l In total, the systematic uncertainty is ~9%, with 
correction factors that range from 1.01 to 1.07. 
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Table 3: Summary of corrections and systematic uncertainties to the ratio of branching 
fractions B{Hb -> HcTr~7T~^Tr~)/B{Hb Hcir~ 



Source 




central value ± syst. error 






W 


B- B° 




Track Reconstruction 




1.00 ±0.04 




IP /Vertex Resolution 




1.02 ± 0.03 




Production/Decay Model 




1.04 ±0.04 




Best Cand. Selection 


1.02 ±0.02 


1.01 ±0.01 1.02 ±0.02 


1.03 ±0.02 


Trigger Efficiency 




1.00 ±0.02 




Fitting 


1.00 ±0.04 


1.00 ±0.04 1.00 ±0.06 


1.00 ±0.04 


Cut on ^Tracks 


0.99 ±0.01 


0.95 ±0.01 0.99 ±0.01 


0.95 ±0.01 


PID 




1.01 ±0.01 




HcD^ background 




0.99 ±0.01 




MC Statistics 


1.00 ±0.04 


1.00 ±0.03 1.00 ±0.04 


1.00 ±0.04 


Total Correction 


1.07 


1.01 1.07 


1.03 


Total Systematic (%) 


8.8 


8.4 10.1 


9.2 



6 Results for Hh HcTT tt+tt 

The results for the ratios of branching ratios are 

B{B^ D+Tx-n+Tx-) 



B{BO D+n 
B{B 



B{B- 
B{Bl Dfn-TT+TT- 
B{B^, ^ Dfn-) 

B{Kl A + TT-TT + TT- 

B{Kl A+TT-) 



2.38 ±0.11 ±0.21 
1.27 ±0.06 ±0.11 
2.01 ±0.37 ±0.20 
1.43 ±0.16 ±0.13, 



where the first uncertainty is statistical and the second is systematic. These measurements 
are all substantially more precise than the current world average values. Naively, one 
might have expected the four branching fraction ratios to be nearly equal. The observed 
differences may be explained in terms of the contributing Feynman diagrams. From Fig. [H 
we see that the primary contribution to B^ — )■ D^'k~{'k^'k~) and -B^ — )■ D+7r~(7r+7r~) is 
from a single decay diagram, an external tree diagram. On the other hand the B~ — )■ 
D^Tx^ [ti'^'k") and A° — >■ A+7r~(7r^"7r~) amplitudes receive contributions from both external 
and color-suppressed tree diagrams. This would suggest that the interference tends to 



be more constructive in i? — >■ D^tt and A? — >■ Aj"7r than in 5 — )■ D^tt~ 



vr ' TT 



and 



Aj 



TT TT TT 



respectively. The role of the various contributing topological amphtudes 
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and the strong phases in i? — > Dtx is discussed in the hterature [12] . In general we see the 
branching fractions for the HcTtttti final states are at least as large or even twice as large 
as the single- TT bachelor states. 

7 Kinematic Distributions and Mass Spectra in the 
7T~7T~^7T~ System 

Since we rely on MC simulation to estimate signal efficiencies, we now compare a few 
distributions between signal MC simulation and data. The higher signal yield — )■ 
D^Ti~ and — )■ D^t[~ii^t[~ decay modes are used, and for each we perform a sideband 
subtraction, where the signal region includes candidates within 50 MeV/c^ of the B^ 
mass, {mso) |15], and the sidebands 60 < |M — m^ol < 110 MeV/c^. For both data and 
simulation, we require events to pass any LO trigger, and signal candidates must satisfy 
the HLTl and HLT2 triggers described in Sec. [H Clearly, two of the most important 
quantities used in our candidate selection are the pt and IP of the daughters from the 

and the recoiling pion(s). Figure H] compares the pt and IP distributions of the 
daughters in data to those from signal MC simulation. Figure |5] shows the corresponding 
comparisons for the recoiling pion(s) in the respective B decay. Overall, the agreement 
between data and MC simulation is very good. 

It is also interesting to examine the 7r~7r"'"7r~ invariant mass spectra for the four signal 
decay modes. Here, we use the sPlot method [I9] to obtain the underlying signal spec- 
tra, based on the event-by-event 6-hadron mass signal and background probabilities. The 
7r~7r'^7r~ mass spectra are shown in Fig. |6l along with signal MC shapes that are normal- 
ized to the same yield as data. We also show several resonant contributions: Di(2420)"'" 
(2%), L)i(2420)° and L)*(2460)° (14% in total), Ac(2595)+ and Ae(2625)+ (9% total), and 
E° and S^"*" (12% total), where the quantities in parentheses are the normalizations rel- 
ative to the total (see Sec. [HI) A prominent structure at low mass, consistent with the 
ai(1260)~ is evident for all decay modes, along with a long tail extending to 3 GeV/c^. 
In all cases, the Sir mass spectrum appears shifted toward lower mass as compared to the 
MC simulation. The simulated value for the ai(1260)~ mass is 1230 MeV/c^, which is 
equal to the central value given in the PDG [TSj of (1230 ±40) MeV/c^. Besides having a 
large uncertainty, the mass as obtained by experiment may be process-dependent, so it is 
difficult to draw any definitive conclusion from this shift. Since both the reconstruction 
and trigger efficiency are flat through this mass region, this small shift in mass does not 
introduce any signiflcant systematic uncertainty in the branching fraction measurement. 

We have also looked at the di-pion invariant masses within the Sir system, shown for 
B^ —7- D^iT~ir~^ir~{a.,h) and -B~ — ?■ D''7r~7r"'"7r~(c,d) in Fig. [71 Contributions from the 
narrow excited charm states, which are discussed in Sec. [HI are excluded. In all cases, in 
the low M{ir~ir^ii') mass region, we see a dominant p^iT~ contribution, consistent with 
the ai(1260)~ resonance. In the higher M{irinT) regions there appears to be an additional 
resonant structure, consistent with the /2(1270) state, in addition to the p° contribution. 
Similar spectra are found for B^ — >■ D'^n^n^n^ and A° — )■ A^7r~7r"^7r~ (not shown.) 
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Figure 4: Comparisons of the and IP spectra for the daughters from the in — )■ 
D^7r~ [(a) and (b)], and from the in — >■ D"'"7r~7r"'"7r~ [(c) and (d)]; Points with 
error bars are data and the sohd hues are simulation. 

The /2(1270) has been previously seen in — )■ D*+7r^7r+7r^ [2T]. The like-sign di-pion 
invariant mass spectra do not show any resonant features. 

8 Contributions from Excited Charm Hadrons 



Within the Hb HcTv^tv^tt 
Ac(2625)+ and S^'^ 



final state, we search for L)i(2420), L'|(2460), Ac(2595)^ 



which may decay to D or A+ with an accompanying tt or tttt 
pair. To search for H* — )■ H^'n'^'n'^ intermediate states, we select events in the 6-hadron 
signal region (±60 MeV/c^ around the nominal mass) and compute the invariant mass 
difference AMtttt = M{HcTi^ti~) — M{Hc) (two combinations per 6-hadron candidate.) 
For the A° ^ S°'++7r±7r-, 2°'++ ^ A+7r±, we use AM^ = M{Hc7r^) - M{Hc) in a 
similar way (one (two) S^"*" (S°) candidates per A° decay.) We also have looked in the 
upper mass sidebands, and the AM^.^- and AM^ distributions are consistent with a smooth 
background shape with no signal component. We look at all data, irrespective of trigger. 
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(GeV/c) IP (mm) 



Figure 5: Comparisons of the px and IP spectra for the bachelor pion in — D^n^ [(a) 
and (b)], and for the 3 pions in — )■ D^ti~tt^'k~ [(c) and (d)]. Points with error bars 
are data and the sohd hnes are simulation. 

to establish signal significances, but for the branching fraction measurement, we use the 
same trigger requirements described in Sec. [71 We choose only one candidate per event 
using the same criteria as discussed previously. We normalize the rates to the respective 
inclusive H^y — )■ Hc7r~7T^7r~ decay, using the same trigger selection as above. We show 
only the AM^^^ and AM^r distributions after the specified trigger, since the distributions 
before the trigger are quite similar, except they typically have 25—30% larger yields than 
the ones shown. 

The AM^TT distributions for and B"^ are shown in Fig. [8] and the AM^r for are 
shown in Fig. [HI For B^, the size of the data sample is insufficient to observe the excited 
Ds states in these hadronic decays. 

Signal yields are determined using unbinned extended maximum likelihood fits. Start- 
ing with 5° (Fig. IHt^a)), we see an excess at AM^^r ~ 560 MeV/c^, consistent with the 
Di(2420)"'". We fit the distribution to the sum of a signal Breit-Wigner shape convoluted 
with a Gaussian resolution, and an exponential background shape. The full width is fixed 
to 25 MeV/c^ [I^ and the mass resolution is set to 7.5 MeV/c^ based on simulation. The 
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Figure 6: Invariant mass of the Svr system in — D"*"?? 7r"''7r (top left), B — )■ 
£)°7r-7r+7r- (top right), S° D+tt-tt+tt- (bottom left) and A+7r"7r+7r" (bot- 

tom right) decays. The data are the points with error bars and the simulations are the 
solid lines and shaded regions. 



fitted yield is 33 ± 8 events and the fitted mean is (562 ± 4) MeV/c^, consistent with 
the expected value. If the width is allowed to float, we find (22.7 ± 8.0(stat)) MeV/c^, 
also in agreement with the world average. Prior to applying the specific trigger selection, 
we find 40 ± 9 signal events, corresponding to a statistical significance of 6.8 standard 
deviations (for one degree of freedom) as determined from the difference in log-likelihoods, 
y/—2ALL, where the difference is taken between the signal yield taken as a free parameter 
and fixed to zero. 

The AMtttt distributions for B~ displayed in Fig. |8]^b) show not only the Di(2420)'', 
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Figure 7: tt+tt invariant mass (2 combinations per candidate) in the Stt system 
for 5° ^ D+TT-TT+TT- when (a) M(7r-7r+7r~) < 1.5 GeV/c^ and (b) M(7r-7r+7r-) > 
1.5 GeV/c^. The corresponding plots for B~ — D^vr^Tr'^Tr" are shown in (c) and (d). 



but also a shoulder at ~600 MeV/c^, consistent with the D2{24Q0)^. Hence, we allow 
for both Di(2420)° and (2460)° signal components, and fix their full widths to the 
PDG values [I5] of 20.4 MeV/c^ and 42.9 MeV/c^ respectively. The means and yields 
are left as free parameters in the fit. The fitted Di(2420)° and 1^2(2460)° yields are 
124 ± 14 and 49 ± 12, with masses that are consistent with the expected values. The 
respective signal yields before the trigger requirement are 165 ± 17 and 63 ± 15 events, 
with corresponding statistical significances of 10.5 and 5.5 standard deviations for the 
Di(2420)° and 1^2(2460)°, respectively. These and B~ decays have also been observed 
by Belle [22]. 

We have also measured the relative fractions of Z}i(2420)° and £'2(2460)'' that do or 
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Figure 8: Invariant mass difference M{D7i' 



candidates, (b) B — )■ D^tt tt^tt signal candidates, (c) B 



M{D), for (a) 5° D+7i-7i+n~ signal 



TT 



through a D* 



intermediate state, and (d) B~ — )■ 7r~-7r"''7r~ not through a D*~^ intermediate state. The 
signal components are the white region (and lightly shaded regions for B~ D^tt'tt^tt^), 
and the background component is the darker shaded region. 



do not decay through D*~^ by taking the subset of candidates with M(D°7r"'") — M{D^) < 
150 MeV/c^ or M(L'°7r+) -M(L'°) > 150 MeV/c^, respectively. The corresponding AM^^ 
distributions are shown in Fig. [Ht^c) and Fig. Et^d). A fit is made to the data as discussed 
previously, and the yields are summarized in Table HI 

For A° (see Fig. [n](a)), we find two well-separated peaks in the AM^t^ distribution, one 
at ~307 MeV/c^, and a second at ~340 MeV/c^, consistent with the expected values for 
the Ac(2595)+ and Ac(2625)+, respectively. The full width of the Ac(2595)+ is fixed to the 
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PDG value of 3.6 MeV/c^, and the mass resolution for each peak is fixed to 2.0 MeV/c^, 
as determined from simulation. The fitted signal yields are 9.7 ± 3.5 and 9.3 ± 3.2 for the 
Ac(2595)"'" and Ac(2625)"'", respectively. Before the trigger, we find signal yields of 10.6 ± 
3.8 for Ac(2595)"'" and 15.7 ± 4.1 for Ac(2625)"'~, corresponding to statistical significances 
of 4.3 and 6.6 standard deviations. Thus we have evidence for A° — >■ Ac(2595)+7r~ and 
observation of A^ — >■ Ac(2625)+7r~. The systematic uncertainties do not change this 
conclusion. These decays have also been reported by CDF |23], but are not yet published. 
The fitted AM^^ values of (306.7 ± 1.1) MeV/c^ and (341.7 ± 0.6) MeV/c^, for the 
Ac(2625)"'" and Ac(2625)+, respectively, are consistent with the known mass differences [15] 
for these excited states. 

We also observe the decays Aj] S°'++7r=F7r", with S° A+vr" or S++ A+7r+. The 
AMtt distributions are shown in Fig.[9]^b-d) for both E° and S^"*" candidates, S° candidates 
only (c), and (d) candidates only. The data are fit to the sum of a Breit-Wigner shape 
convolved with a Gaussian resolution function and a smooth threshold function. The full 
width is fixed to 2.2 MeV/c^ [13] in all cases, and the AM^r resolution is fixed to 1 MeV/c^ 
based on simulation. The combined S° and signal has a statistical significance of 6.0 
standard deviations. The S° and S^"*" signals have statistical significances of 4.9 and 3.5, 
respectively. These decays have also been seen by CDF [23] . 

Table S) summarizes the yields for the various excited charm states for both the full 
data sample and after the trigger selection as well as the yields in the normalizing modes 
(after trigger selection.) 

Table 4: Summary of yields for the signal and normalization modes. Below Di and D2 
refer to the Di(2420) and D|(2460) mesons, respectively. 



Decay 


H*7t{7t) Sk 


j;nal Yields 


HcTl TT+TT 




All 


Trig. Sel 


Trig. Sel 


S° ^ Dfir-, Dt D+TT-TT+ 


41 ±8 


33 ±7 


1741 ±55 


B- D°7r-, D°TT-TT+ 


165 ± 17 


126 ± 14 


1386 ±51 


B- D^TT-, ^ D*+TT- 


111 ± 14 


75± 12 


1386 ±51 


B~ D^^TX^, D° D^7r~7i+, non-D* 


57± 10 


52 ±9 


1386 ±51 


B- Dfn~, Df D^n-n+ 


66± 15 


49± 12 


1386 ±51 


B~ DfTT-, Df D*+Tx- 


46± 12 


34± 10 


1386 ±51 


B~ DfTT-, Df D^n-7r+, non-D* 


23 ±9 


18 ±8 


1386 ±51 


AO ^ Ac(2595)+7r- 


10.6 ±3.8 


9.7 ±3.5 


312 ±23 


AO ^ Ac(2625)+7r- 


15.7 ±4.1 


9.3 ±3.2 


312 ±23 


AO ^ S0'++7rT7r- 


29.3 ±7.0 


24.9 ±6.2 


312 ±23 


AO ^ S07r-7r+ 


19.6 ±5.7 


16.2 ±5.0 


312 ±23 


A" S++7r-7r- 


10.1 ±4.0 


9.3 ±3.7 


312 ±23 
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Figure 9: Intermediate resonances contributing to the A° — )■ A+tt TT'^ir decay. Shown are 
distributions for (a) M(A+7r-7r+) - M(A+), with Ac(2595)+ and Ac(2625)+ contributions, 
(b) M(A+7r^) - M(A+) (3 combinations per A° candidate) (c) M(A+7r-) - M(A+) (2 
combinations per A° candidate), and (d) M(A+7r+) — M(A+) (1 combination per A^ 
candidate), showing the intermediate Sc states. The hne is a fit as described in the text, 
and the shaded region is the fitted background. 



The branching ratios for these modes are computed using: 

B{H, ^ H^nin)) x B{H: ^ HM^)) _ , , , . -i 

where if* refers to one of the observed excited charm states, A^'signai and A'norm are the 
number of reconstructed decays in the signal and normahzation modes after the trigger 
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requirement, 6^°} is the reconstruction and selection efficiency relative to the normalization 
mode, and e];^|g|gg[ is the relative trigger efficiency. All efficiencies are given for the mass 
region 0.8 GeV/c^ < M{tt-tt+tt-) < 3 GeV/c^. 

The relative reconstruction, selection and trigger efficiencies, shown in Table [5|, are 
evaluated using MC simulations. The Z}i(2420)0 and 1)2(2460)° are each assumed to 
decay 70% through D*+Tr- D^tt+tt- and 30% non-resonant D^tt+tt-. The Di{2A20)+ 
is taken to be 100% non-resonant D~^n~7t'~^. The Ac(2595)"'" decay is simulated as 36% 
36% S+^TT" and 28% non-resonant A+tt"?!"*". The Ac(2625)"'" decay is assumed 
to be 100% non-resonant A^7r^7r+. The Sc(2544) baryons are simulated non-resonant in 
phase space. 

The relative efficiencies agree qualitatively with our expectations based on the kine- 
matics and proximity to threshold for these excited charm states. The differences in the 
relative efficiency between the pairs of excited charm states for a given 6-hadron species 
are negligible compared to the uncertainty from our limited MC event sample, and we 
use the average relative efficiency for each pair of decays. 

Table 5: Summary of the relative reconstruction and selection efficiencies (e^^}) and trigger 
efficiencies (e^^|g|ggi) for the excited charm hadron intermediate states with respect to the 
inclusive HcH'n^ii' final states. Below Di and refer to Di(2420) and 1^2(2460), 
respectively. The uncertainties shown are statistical only. 



Decay 




_rcl 
^scl 

(%) 


trig|sel 
(%) 


_rcl 
^total 

(%) 






0.83 ±0.06 


1.05 ±0.09 


0.87 ±0.10 


B- -> {Dl Df)n- 




0.70 ±0.04 


1.24 ±0.07 


0.86 ±0.07 


B- (DO, D*0)7r" 


(via D*) 


0.66 ±0.05 


1.29 ±0.08 


0.84 ±0.08 


B- ^ (DO, D*0)vr- 


(non-D*) 


0.78 ±0.06 


1.15 ±0.10 


0.91 ±0.11 


AO ^ (A,(2595), Ae(2625)+)7r- 


0.52 ±0.03 


1.30 ±0.07 


0.67 ±0.06 


AO ^ S^+TTTT, SO--^ 


'+ ^ A+vr^ 


0.67 ±0.05 


1.10 ±0.13 


0.75 ±0.10 



The dominant sources of systematic uncertainty are the limited MC sample sizes 
and the ffi model. Starting with the B^, the uncertainty due to limited MC statis- 
tics is 11%. For the ffi model, the largest source of uncertainty is from a possible 
D2(2460)+7r~, ^2(2460)+ -> D+7r"7r+ contribution. If this contribution is included in 
the ffi using a Breit-Wigner shape with mean and width taken from the PDG [15], the 
returned signal yield is O^q. If we assume isospin symmetry, and constrain this fraction 
(relative to Di(2420)) to be (40 ±11)%, the ratio found for the B' decay, the ffited 
fiO Di(2420)+7r-, Di(2420)+ D+Tr~n+ signal yield is 26 ± 6 events. We take this 
as a one-sided uncertainty of ^Sw- Sensitivity to the background shape is estimated 
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by using a second order polynomial for the background (3%). The mass sidebands, 
which have a i5i(2420)+ fitted yield of 2^2 events from which we conservatively assign 
as a one-sided systematic uncertainty of For the signal decays, 4% of events have 
M(7r~7r+7r^) > 3 GeV/c^, whereas for the -Di(2420) + , we find a negligible fraction fail 
this requirement. We therefore apply a correction of 0.96 ± 0.02, where we have taken 
50% uncertainty on the correction as the systematic error. The systematic uncertainty 
on the yield in the — )■ Z}+7r~7r+7r~ normalizing mode is 3%. We thus arrive at a total 
systematic error on the B'^ branching fraction ratio of 

For the B~ , we have a similar set of uncertainties. They are as follows: MC sample 
size (8%), background model (1%, 2%), L)i(2420)° width (2%, 4%), L)*(2460)° width 
(1%, 3%), where the two uncertainties are for the (Di(2420)°, D2(2460)°) intermediate 
states. We have not accounted for interference, and have assumed it is negligible compared 
to other uncertainties. A factor of 0.98 ± 0.01 is applied to correct for the fraction 
of events with M(7r~7r+7r^) > 3 GeV/c^. Including a 3% uncertainty on the B^ — j- 
D^TT~7r^7T~ yield, we find total systematic errors of 9% and 10% for the Di(2420)'' and 
1^2(2460)'' intermediate states, respectively. For the D* sub-decays, the total systematic 
uncertainties are 10% and 11% for 5" ^ Di(2420)°7r-, Di(2420)° ^ D*+7r- and B~ 
L'2(2460)°7r~, D;(2460)° D*+7r~, respectively. For final states not through D*, we 
find a total systematic uncertainty of 13% for both intermediate states. In all cases, the 
dominant systematic uncertainty is the limited number of MC events. 

For the A° branching fraction ratios, we attribute uncertainty to limited MC sample 
sizes (8%), the A+(2595) width AO ^ A+tt-tt+tt- signal yield (3%), and apply a 

correction of 0.96 ± 0.02 for the ratio of yields with M(7r~7r+7r~) > 3 GeV/c^. In total, 
the systematic uncertainties on the A+(2595)"'" and Ac(2625)"*" partial branching fractions 
are 'tHy^ and iblO%, respectively. 

For the Sc'"*"^ intermediate states, the systematic uncertainties include 14% from finite 
MC statistics, and 4% from the width. For the Sc'^"*" simulation, 10% of decays 

have M(7r~7r+7r~) > 3 GeV/c^, compared to 4% for the normalizing mode. We therefore 
apply a correction of 1.06±0.03 to the ratio of branching fractions. All other uncertainties 
are negligible in comparison. We thus arrive at a total systematic uncertainty of 16%. 
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The final partial branching fractions are 



B{B- D?7r+, ^ D^-n^ 



B- DH-n+Tx- 
B{B- DIti+, D\ D*+7i- 



B- DOtt'tt+tt- 



B- DH-Ti+Ti- 
B{B- -> D^°7r+, D*^ D^Tv-TT^ 



B- DH-TT+TT- 

B{B- 0^71+, Df D* 



TV 



B- ^ DH-TT+TT^ 

B{B- ^ D*07r+, Df ^ D07r-7r+)„o„_z) 



B- D^TT-n+TT- 



B{kl Ac(2595)+7r+, Ac(2595)+ ^ A+7r-7r+ 



i3(A0 -> Ae(2625)+7r+, Ac(2625)+ ^ A+7r-7r^ 



A^ A + TT-TT + TT- 



TT ' TT 



^ A: 



TT 



A[J A + TT-TT + TT- 



B{At ^ E++7r-7r-, S++ ^ A+ 



A^ A+TT-TT + TT- 



< 



2.1±0.5^°J)% 
10.3 ± 1.5 ±0. 

9.3 ± 1.6 ±0.9)% 
4.0 ±0.7 ±0.5)% 
4.0 ± 1.0 ±0. 
3.9± 1.2±0. 

1.4 ±0.6 ±0.2)% 
!.0% at 90% C.L. 



4.4 ± 1.7+^ 



OAJ 



4.3 ± 1.5 ±0.4)% 
11.4 ±3.1 ± 1.8)9^ 

7.4 ± 2.4 ± 1.2)% 
4.2 ± 1.8 ±0.7)%, 



where the first uncertainties are statistical and the second are systematic. For the modes 
with D*+, we include a factor B{D*+ D^n+) = (0.677 ± 0.005) [15] to account for 
unobserved D*~^ decays. The first four and the sixth of these decays have been pre- 
viously measured by Belle [22] with comparable precision. To compare our results to 
those absolute branching fractions, we multiply them by the relative B^ [B~] branch- 
ing fractions in Eq. [H and then in turn by — )■ £)+7r") = (2.68 ± 0.13) x 10~^ 
[B{B^ —7- D^Ti^) = (4.84 ± 0.15) x 10^^.] The resulting absolute branching fractions are 
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^ Di(2420)-7r+, L)i(2420)- ^ D+7r-7r+) = (1.3 ± 0.3+°;^) x 10"^ 
B{B- L)i(2420)°7r+, L)i(2420)° ^ L)°7r-7r+) = (6.3 ± 0.9 ± 0.9) x 10"^ 
B{B- Di(2420)°7r+, Di{2A20f D*+7r") = (5.8 ± 1.0 ± 0.9) x 10"^ 
B{B- L)i(2420)°7r+, L)i(2420)° ^ D\+7r-)non~D' = (2.5 ± 0.4 ± 0.4) x 10"^ 
B{B- ^ D;(2460)°7r+, D;(2460)° ^ D*+7r") = (2.5 ± 0.7 ± 0.4) x 10"^ 



where the uncertainties are statistical and total systematic, respectively. The correspond- 
ing values obtained by Belle are: (0.89;[|;|) x 10-^ (6.5+};^) x 10-^ (6.8 ± 1.5) x 10-^ 
(1.9+o;6) X 10"^, and (1.8 ± 0.5) x 10"^ 122]. Our resuhs are consistent with, and of 
comparable precision to, those measurements. 

Preliminary results on the A+(2595)+7r-, A+(2625)+7r- and A|^ 

S^'+^TT^vr" decays have been reported by CDF [23]. Our values are consistent with 
these (unpublished) results. 



9 Summary 

In summary, we have measured the branching fractions for Hb — )■ Hc'k~ti'^ti~ decays 
relative to Hh — )■ HcT^~ . The ratio of branching fractions are measured to be 



BiB^ D+TT-TT+TT-) 

^ ^ = 2.38 ±0.11 ±0.21 



B{B^ D+7T-) 
B{B- D^n-TT+TT- 
B{B- D^Ti-) 

B{B^^ D+TT-TT+TT- 

i3(50 ^ D+VT") 
g(A° ^ A + TT-TT + TT- 

B{Kl ^ A+TT-) 



: 1.27 ±0.06 ±0.11 
: 2.01 ±0.37 ±0.20 
1.43 ±0.16 ±0.13. 



At low 37r mass, these decays appear to be dominated by the ai(1260) resonance. We 
have also measured several partial decay rates through excited charm states. The yields of 
Hiy — 7- iJc^ TT+TT" relative to — >■ Hc7r~ are in the range of 20—40%. If the relative rates 
in the Cabibbo-suppressed decays, such as B^ — )■ DfR^n^ii^ and B^ — t- DK~ti^ti~ 
relative to 5° — > DfK^ and B" — )• DK~, respectively, are comparable, they could be 
useful for measuring the weak phase 7. 
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